Plaquette bond order wave in the quarter-filled extended Hubbard model on the 

checkerboard lattice 
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An extended Hubbard model (including nearest-neighbor repulsion and antiferromagnetic spin 
exchange) is investigated on the frustrated checkerboard lattice, a two-dimensional analog of the 
pyrochlore lattice. Combining Gutzwiller renormalized mean-field (MF) calculations, exact diago- 
nalization (ED) techniques, and a weak-coupling renormalization group (RG) analysis we provide 
strong evidence for a crystalline valence bond plaquette phase at quarter-filling. The ground state 
is twofold degenerate and breaks translation symmetry. The bond energies show a staggering while 
the charge distribution remains uniform. 
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I. INTRODUCTION 

Many frustrated quantum magnets exhibit complex 
behavior deeply rooted in the macroscopic degeneracy of 
their classical analog. Spin liquids characterized by the 
absence of any kind of symmetry breaking and valence 
bond crystals (VBC) characterized by translation sym- 
metry breaking only (SU(2) symmetry is preserved) are 
exotic phases possibly realized in some two-dimensional 
(2D) or three-dimensional (3D) frustrated magnets. 1 The 
VBC instability in 2D and 3D are higher dimensional re- 
alizations of the one-dimensional (ID) spin-Peierls mech- 
anism which is known to occur, e.g., in the Heisen- 
berg spin-1/2 chain with nearest and next-nearest neigh- 
bor exchange interaction above some critical frustration 
(even in the absence of lattice coupling). 2 In that case, 
the chain spontaneously dimerizes leading to a twofold 
degenerate gapped ground state (GS). Similarly, the GS 
of the Heisenberg model on the checkerboard lattice (the 
square lattice with diagonal bonds on every second pla- 
quette, see Fig. [1]) was shown to be twofold degenerate 
with long range plaquette VBC order characterized by a 
Q = (tt, it) ordering wavevector£ In this state the spins 
pair up in four-spin singlet units located on every sec- 
ond void plaquette giving rise to stronger bonds on those 
plaquettes. 

So far the investigation of charge degrees of freedom 
in geometrically frustrated systems has been poorly ex- 
plored. Most computations deal with a small fraction 
of holes doped into a Mott insulator— or with itinerant 
systems at half-filling. 5 

It has been suggested recently^ that for strongly cor- 
related electrons on frustrated lattices such as the 2D 
kagome or the 3D pyrochlore lattice a related sponta- 
neous symmetry breaking away from half-filling might 
occur at the fractional filling n — 2/3 and n — 1/2, re- 
spectively. These so-called bisimplex lattices are arrays 
of corner sharing triangular or tetrahedral units. Inter- 



estingly, although the GS is twofold degenerate due to the 
breaking of inversion symmetry that leads to two types 
of non-equivalent units, lattice translational symmetry is 
preserved and all sites remain equivalent, i.e., there is no 
charge density wave (CDW) ordering. Each site of the 
lattice belongs to two non-equivalent units and the two 
types of units have different bond strength, i.e., different 
expectation values for the kinetic and magnetic energy, 
which is the typical characteristic of a bond order wave 
(BOW). 

The checkerboard lattice can be viewed as a 2D corner- 
sharing array of tetrahedra. As shown in Ref. the t-J 
model with two electrons on a single tetrahedron has a 
very robust non-degenerate and non-frustrated ground 
state which is given by an equal amplitude superposition 
of the six states having a singlet on one bond. Therefore, 
it was conjectured in Ref. [6| that a BOW might occur 
quite generally on a bipartite and corner-sharing arrange- 
ment of tetrahedra at quarter-filling, i.e., that the BOW 
instability occurs not only in the 3D pyrochlore lattice 
but also in the 2D checkerboard lattice. 

In the present paper, we provide numerical and ana- 
lytical evidence for this new type of spontaneous BOW 
instability for correlated electrons on the checkerboard 
lattice at quarter-filling [n = 1/2). The spontaneous 
symmetry breaking produces a doubling of the unit cell 
like in the spin-Peierls and VBC scenarios and does not 
lead to a breaking of inversion symmetry as on the py- 
rochlore lattice. The BOW instability occurs on the 
checkerboard lattice at quarter-filling in the presence of 
nearest-neighbor (n.n.) interactions (repulsion and mag- 
netic exchange). 

The BOW phase bears similarities with the VBC phase 
at half- filling (Heisenberg model), namely long range 
plaquette-plaquette correlations with the same wavevec- 
tor, but the spin singlets are formed on the crossed pla- 
quettes for the BOW at quarter-filling whereas they are 
formed on the empty plaquettes for VBC at half-filling. 
We report the results of strong-coupling (U = oo) ap- 
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proaches in Sec. [TTJ where Sec. Ill Al is dedicated to a 
Gutzwiller MF calculation and Sec. Ill Bl contains ED re- 
sults. 

In Sec. IIIII we derive an effective weak-coupling model 
on the square lattice and analyze it with RG methods. 
The Sec. IIII Al contains a MF and two-patch RG analysis 
and in Sec. IIII Bl the results of the more sophisticated 
A-patch RG are presented. 

The different approaches provide a very consistent pic- 
ture and show that the BOW instability on the checker- 
board lattice at quarter-filling occurs both for weak and 
strong coupling and for quite general repulsive interac- 
tions. 
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A. The Hamiltonian 

We study the Hamiltonian H = Ho + H lnt at quarter- 
filling (n = 1/2) on the checkerboard lattice. The kinetic 
part is given by 

#o = -*E£ (4^+ h - c -) a) 

(ij) °-=U 

with the positive hopping matrix element t. The sum 
Y^tUj) runs over au bonds on the checkerboard lattice. 
The interaction part of the Hamiltonian is given by 

H^t = U/, n^riji + J}^ Sj ■ S j + V mnj (2) 

» (ij) (ij) 

with onsite repulsion U, n.n. repulsion V, and n.n. spin 
exchange J. For U = oo we obtain the strong coupling 
Hamiltonian 

H t .j_ v = VHoV + J S i ■ S j + V n * n i ' ( 3 ) 

(ij) (ij) 

= Ht.j + V'Y.mnj, (4) 

(ij) 

where V is the projection operator that enforces the sin- 
gle occupancy constraint and V' = V + J/4. For V' = 
the strong coupling Hamiltonian reduces to the usual t-J 
model. 



II. STRONG COUPLING 

A. Gutzwiller renormalized mean-field 

We study the t-J-V model Hamiltonian on the 
checkerboard lattice. This Hamiltonian describes the 
physical processes in the restricted Hilbert space of con- 
figurations with no doubly occupied sites. Physically, we 
expect V > J but we do not necessarily assume V large 
compared to t. The local constraints of no doubly occu- 
pied sites are replaced by statistical Gutzwiller weights 7 



FIG. 1: (Color on-line) Schematic pattern of the plaquette 
phase. Different line styles (colors) correspond to different 
electron hopping amplitudes and spin-spin correlations (see 
Fig.©. 

and a decoupling in the particle-hole channel leads to the 
following MF Hamiltonian, 

H M F= - tJ2J29ij( C ia C j* + h-C.) 

(ij) a 

Yu (l4 J + v ) (xji c l c j* + h - c - - \x» i 2 ) 

(ij)o- 

+ Vj2(( n t) n i + (nj)fH - (ni){nj)), (5) 
(ij) 

where the Gutzwiller weights have been expressed in 
terms of local fugacities Zi = (1 — (rij))/(l — (rii)/2), 
gjj = ^/ZiZj and gfj = (2 — 2j)(2 — Zj), to account for pos- 
sible (small) non-uniform charge modulations (if any). 8 
The self-consistency conditions are implemented as Xji = 
(cj a c ia ). Recently, this approach was successfully used to 
investigate the properties of 4 x 4 checkerboard-like BOW 
in the lightly doped t-J model on the (non- frustrated) 
square lattice.— 

The MF equations are solved self-consistent ly on 16 x 
16 and 48 x 48 lattices starting from a random bond and 
site configuration. We assume a fixed physical value of 
J/t = 1/3 and systematically vary the ratio V/t. We 
find that in the range 0.3 < V/t < 1.75 the MF equa- 
tions converge to a unique solution given by the checker- 
board pattern shown in Fig. [T] Interestingly, finite size 
effects are weak for this range of parameters. As far as 
symmetry properties are concerned, the modulated struc- 
ture found here corresponds to a doubling of the unit cell 
with two non-equivalent "weak" and "strong" tetrahedra 
(crossed plaquettes) hence leading to four types of bonds 
schematically represented by different line styles (colors) . 
Each bond is characterized by the expectation values of 
the hopping and of the spin-spin correlations which are 
plotted in Fig. [5] It is important to emphasize here that 
no particular superstructure was assumed in the calcula- 
tion (the equations are solved in real space starting from 




0.15- 



0.10 



(a) N = 20 




0.05*-^ 



o 



CD 



O 
O 



Q_ 

o 



0.5 1 1.5 2 0.5 1 1.5 2 
V/t V/t 

FIG. 2: (Color on-line) Expectation values vs. V/i of the hop- 
ping \{(cl a c jcr +h.c.))| (left) and exchange |(S, ■ Sj)\ (right) 
operators on the four non-equivalent bonds of the plaquette 
phase (see Fig. [T|. A different line style and color is used 
for the different plaquettes and the thicker (thiner) lines cor- 
respond to the outer (crossed) bonds of the plaquette. The 
results are obtained from an unrestricted MF calculation on 
a 48 x 48 lattice for J/t = 1/3. 



a random state) and that our solution space a priori al- 
lows for CDW's and/or bond currents (complex values 
for Xij)- However, our solution exhibits real values of 
Xij (it does not break time-reversal symmetry) and has 
a uniform charge density on all sites. It is therefore a 
pure BOW. Note that two types of bonds (with slightly 
different amplitudes) are found within each tetrahedron. 
Although the MF treatment, very likely, overestimates 
the amplitude of the modulations^ these results con- 
vincingly suggest a quite strong BOW instability of the 
checkerboard lattice. Note also that for large V/t val- 
ues our MF approach fails to converge as the system of 
non-linear equations becomes unstable and oscillates be- 
tween two attractors. In fact, in the strong V/t limit, 
the "tetrahedron rule"— (or ice rule) should be obeyed, 
namely there should be exactly two particles within each 
crossed plaquette. It is clear that the simple decoupling 
scheme of the V term does not fulfill this "hard-core" 
constraint. 



B. Exact Diagonalizations 

We study the strong coupling t-J-V model ^ within 
ED on finite-size checkerboard samples of N = 16, 20, 
and 24 sites at the filling n = 1/2 J£ 

In finite, periodic systems the BOW instability can 
be detected using a correlation function of the bond 
strengths. Here we chose to work with the kinetic term. 
The correlation function is defined as: 



C K [(i,j),(k,l)] = (lC(i,j) AC(fc,Z)> 



(JC(i,j))(JC(k,l)), 
(6) 
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FIG. 3: (Color online) Correlation function of the kinetic en- 
ergy (Eq. [6} of (a) a 20 sites and (b) a 24 sites checkerboard 
sample at n = 1/2, J/t = 1/3, and V' /t = 7. The black, 
empty bonds denote the same reference bond, the red, full 
bonds negative and the blue, dashed bonds positive correla- 
tions. The line strength is proportional to the magnitude of 
the correlations. 



where 



h.c. 



(7) 



and is a n.n. bond while (k, I) denotes any other 

bond which is linked in the Hamiltonian. We evaluate the 
correlation function for all possible combinations without 
common sites in the ground state of several samples at 
n= 1/2 and J/t = 1/3, for varying V'/tM Note that a 
BOW of the type presented in Section Til Al is character- 
ized by long range correlations of C/c corresponding, in 
the thermodynamic limit, to a modulation of (JC(i.j)) to 
be compared directly to the MF values in the left panel 
of Fig. [2] (the long distance correlation function is pro- 
portional to the magnitude-squared of the kinetic energy 
modulation). 

First we plot a real space picture of the correlations on 
the N = 20 and N = 24 samples at a value of V/t = 7. 
As discussed below this is the value of V/t where the 
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FIG. 4: Behavior of the absolute value of a particular bond 
correlation function (Eq.|6} at largest distance as a function of 
V'/t at fixed J/t = 1/3 and n = 1/2 for several system sizes. 
Correlations are weak at small V'/t, peak around V'/t ~ 7 
and then decrease again towards the extremely large V'/t 
limit, where another phase might appear. 



correlations are strongest. The reference bond uniquely 
belongs to a certain class of crossed plaquettes. Based on 
the theoretical picture one expects the correlation func- 
tion to be positive for all bonds on the same type of 
crossed plaquettes and negative on the others. This is 
indeed what is seen on both samples shown in Fig. [3] 
Note also that the correlations are rather regular and 
uniform throughout the system. 

We now track the evolution of the correlations as a 
function oiV'/t by choosing the value of the correlation 
function at the largest possible distance on a given sam- 
ple. We compare the evolution on the three finite size 
samples in Fig. |4j The correlations all start at small val- 
ues for V'/t < 2, and in that region an inspection of the 
spatial structure similar to Fig. [3] reveals that the bond 
correlation pattern has defects, i.e., correlations with the 
wrong sign. This makes it difficult to conclude firmly on 
the presence of the BOW phase at small repulsion based 
on the ED simulations only. The abrupt jump seen on 
the N = 16 sample is due to a level crossing, which is ab- 
sent on the larger samples. For intermediate values of the 
repulsion the pattern becomes correct on all bonds and 
the amplitudes peak around a repulsion of V^axA ~ ? f° r 
the largest two samples. Beyond this value the correla- 
tions weaken again, but the pattern remains qualitatively 
intact. The physics at very large V/t, where states fulfill- 
ing the "tetrahedron rule" are dominant, is of particular 
interest but beyond the scope of this work. Calculations 
for spinless fermions (i.e. polarized electrons)^ reveal a 
more complicated translation symmetry breaking state 
with a larger supercell. In the case of bosons however, a 
similar plaquette BOW is found although involving the 
void plaquettes where cyclic exchange terms act in this 
limitiS. Preliminary numerical calculations and analyt- 



ical considerations^ show further that this second type 
of BOW phase (not shown) is also realized for fermions 
when spin degrees of freedom are included. The investi- 
gation of the transition occuring between these two types 
of plaquette phases is left for future work. 



III. WEAK COUPLING 

A. Two-patch RG and MF analysis 

We turn now to the weak-coupling analysis of the 
Hamiltonian H = Hq + H mi with Hq and H nrt given 
by Eq. ([T]) and Eq. ©, respectively. The quadratic 
Hamiltonian, Hq, in reciprocal space consists of a dis- 
persing band and a flat band lying on top (for t > 0) 
of the dispersing band. If we introduce a chemical po- 
tential term with /i = —It we can write Hq in the form 

Ek ff Ck7k CT 7 kCT + 4i ^flat where the operator N nat counts 
the number of electrons in the flat band at At. In the 
weak-coupling limit these states do not affect the low 
energy physics of the system and therefore they will be 
dropped in the following analysis. The dispersion of the 
other band is given by £k = — 2t(cos/ci + cos £2) an d is 
identical to the tight-binding dispersion on the square lat- 
tice. It is particle-hole symmetric, perfectly nested with 
the nesting vector Q = (tt,tt), and has a logarithmically 
divergent density of states at the Fermi energy. 

For small couplings U, V, and J we obtain an effective 
interaction Hamiltonian (Appendix [S]) given by 



i a 2 <r ' ^k3 <r ' ' k4 



(8) 



ki...k 4 



with 

3ki...k 4 



v=\ 
2 

E( 4 ^k lk Xk 3 - zjfUJUJ (9) 



with 



= cos(fe„/2), / k V = e£_ k '«, V = V- J/4, 
q = ki + k2 — k3 — k-i, and = q— 2(kj — kj). The two 
saddle points Pi — (0,7r) and P 2 — (n, 0) of the disper- 
sion lead to the logarithmic divergence of the density of 
states. Therefore, we can characterize very weak interac- 
tions by the values of the function .gki...k 4 where all four 
momenta ki • • • k4 lie on one of the two saddle points. 
Adopting the notation of Ref. [TH we have the following 
four coupling constants: 



ffki...k 4 



-2 J ki,k 3 E Pi 

+W k!,k 4 ePi 

ki,k 2 ePi 

U + 2V-J ki . . . k 4 £ Pi 



(10) 



From the RG equations of Ref. [l5| we see that for posi- 
tive values of £/, V, and J the coupling g\ flows to — oo 



and the coupling g 2 flows to +00 j2! whereas the other 
couplings flow to 0. This shows that the charge density 
wave susceptibility is diverging most rapidly under the 
RG flow. As shown in Fig. [5] the CDW instability on the 
square lattice corresponds to the BOW instability on the 
checkerboard lattice, i.e., both instabilities break exactly 
the same symmetries. Note, that within the framework 
of this RG scheme, we can not determine whether a usual 
CDW or a so-called charge flux phase, which is a CDW 
with a d-wave form factor, is stabilized. In the follow- 
ing we will denote these two phases with s-CDW and 
d-CDW. (Note, that for J = only the coupling g 2 di- 
verges. In this case we can not even determine whether a 
CDW of SDW phase is stabilized.) In order to show, that 
the s-CDW phase is in fact favored over the d-CDW, at 
least in a mean-field analysis, we restrict the interaction 
Hamiltonian J5J to the CDW channel and obtain: 



#cdw = Jj^Yl Vkk ' H 7L7k+Q, CT 7k'+Q, CT '^'a' (11) 
with 



kk' 



Vkk' — 2gtk,k'+Q,k',k+Q — ffk,k'+Q,k+Q,k' (12) 

+ (13) 



Vr 



kk' 



Vkk' — — Vo [(1 — cos fei cos fe)(l — cos fci cosfe) (14) 

+ (sin 2 ki + sin 2 k 2 )(sin 2 k[ + sin 2 k' 2 )/2\ 

Vfcfc/ = — Vo(cos k\ — cos /s2)(cos k[ — cos k' 2 ) (15) 

Vo sin ki sin k 2 sin k[ sin k' 2 (16) 



V& = 



u 



V^ k , = + — (sin fci sin k[ + sin k 2 sin k' 2 ) 



(17) 



where Vq = (V+J)/2. In (| 1 5() we dropped a term propor- 
tional to cos ki + cos k 2 , as it vanishes along the FS, and 
in (fT7|) we dropped terms proportional to V/U or J/U. 

With A k = i Ek' VwFu and F k = EM+q^) 
we obtain the linearized gap equation 



tanh(g k V2T c ) 
2&< 



Ak'. 



(18) 



Note, that the d'- (fTflj) and the p-wave (IT71) instability do 
not open a gap at the saddle points, furthermore, the p- 
wave instability is strongly repulsive. For the d-CDW 
state the pair potential (fT5)) is separable and the lin- 
earized gap equation can be written in the simpler form 



1 



^^(COS/C! 
k 



cos k- 



2 tanh(g k /2T c d 
' 2£k 



(19) 



The pair potential for the s-CDW state (|T4|) is not sep- 
arable. But if we neglect for a moment the second line 
in (fT4|) we obtain an analogous expression to (fT9|) for the 
critical temperature of the s-CDW with cos fci — cos k 2 
replaced by 1 — cos k\ cos k 2 . As on the FS we have 
cosfc2 = —cos fci and as 12 cos fci I < 1 + cos 2 fci we 
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FIG. 5: (Color on-line) Correspondence between the s-CDW 
phase on the effective square lattice (left) and the plaque- 
tte phase or BOW phase on the original checkerboard lattice 
(right). 



know that the critical temperature of the s-CDW phase 
is higher than the critical temperature of the d-CDW 
phase. Including the non-negative second term in (|14p 
would only lead to a further increase of the critical tem- 
perature. Note, that the d- and the s-CDW potentials are 
identical for the momenta on the saddle points. There- 
fore it is not surprising that they cannot be distinguished 
by the two-patch RG method. However, in the mean-field 
we can see that the s-CDW is favored over the d-CDW 
state, as it opens a finite gap along the entire FS. It is 
possible to perform such a mean-field analysis also for 
SDW and superconducting instabilities. For supercon- 
ductivity we have strong repulsion (oc U) in the d- and 
in the s-wave channel, and weaker (oc V, J) and mainly 
repulsive interactions in the d'- an the p-wave channel. 
It is therefore quite clear, that the superconducting in- 
stabilities can not compete with the s-CDW instability. 
For the pair potential of the SDW instabilities, we should 
replace only V by V/2 and U by -U in Eqs. !|14H17p . 
For J > we have V = (V + J)/2 > V/2 and therefore 
the s-, d- and d'-CDW instabilities are favored over the 
corresponding SDW instabilities. The p-SDW instability 
is strongly attractive but has the handicap, that it does 
not open a gap at the saddle points, therefore, for weak 
interactions the s-CDW state will still be favored over 
the p-SDW state, i.e., for interaction parameters (U, V, J) 
given by (era, av, aj) we can find for all positive values 
of (u,v,j) an ao such that for < a < ocq the s-CDW 
state is stabilized. 

In conclusion, in this section we showed with RG and 
mean-field arguments, that for weak enough repulsive 
and antiferromagnetic interactions the s-CDW phase, 
which corresponds to a bond order wave state on the 
checkerboard lattice, is stabilized. 



B. iV-patch functional RG method 

As the two-patch model in the last section could not 
determine whether the s- or the d-CDW ground state is 
favored, we advance to a more refined A^-patch functional 
renormalization group (fRG) study of the effective model 
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for the dispersive band. This allows us to study the whole 
Fermi surface and, in particular, the competition between 
s-CDW and d-CDW state. The effective interactions in 
these channels and also the bare effective interaction pos- 
sess a non-trivial k-dependence outside the saddle point 
regions which can be taken into account in the ./V-patch 
fRG scheme. 

We set up a one-loop TV-patch fRG calculation of 
the one-particle irreducible interaction vertex using the 
methods described in Refs. [HI- This fRG scheme inte- 
grates out all intermediate one-loop processes (particle- 
hole diagrams and particle-particle contributions) with a 
decreasing energy scale A. In the one-loop diagrams at 
scale A, one intermediate particle is in the energy shell 
with band energy ±A, while the other particle is further 
away from the Fermi surface. The RG flow is started 
at Ao ~ bandwidth. The integration of the one-loop 
processes in the flow down to lower A renormalizes the 
wavevector-dependence of the effective interaction. From 
the growth of certain components toward low scales, the 
dominant correlations can be read off. In parallel with 
the running interactions, the flow of various static sus- 
ceptibilities can be obtained. These quantities are again 
renormalized by the scale-dependent interactions via one- 
loop processes. Their growth can be compared as a func- 
tion of the RG scale. In the case of a divergence of the in- 
teractions at a nonzero critical scale A c , the fasted grow- 
ing susceptibility determines the dominant instability or 
ordering tendency of the system. 

For the spin-rotationally invariant system, the scale- 
dependent interaction vertex can be expressed using a 
running coupling function VX(k, k',k + q). The initial 
coupling function, V\ (k, k', k + q), is obtained from 
Eq. §§§ as ,9k+q,k'-q,k',k- Here the frequency dependence 
has been neglected as in most previous applications of the 
method for correlated fermions. 17 In order to describe 
the wavevector dependence, we use 32 angular patches 
around the Fermi surface-^ These patches are again split 
up in 3 patches, one above the Fermi surface down to 
band energies of 0.2t, one including the Fermi surface in 
the band energy window between —0.2t and 0.2t, and one 
below this energy window. We have checked our results 
using other partitions. 

We have run the fRG for half band filling and various 
choices of the interaction parameters U, V and J in (jHJ. 
From a perturbative point of view, the model with pure 
onsite repulsion and V = J = is remarkably stable. For 
instance at U = 2t and half filling of the dispersive band, 
the fRG flow of the interactions remains finite down to 
low scales and temperatures ~ 10~ 4 t without any indi- 
cation for an instability. This stability is quite surprising 
for a perfectly nested Fermi surface with van Hove sin- 
gularities at the Fermi level. However, as already found 
in the two-patch model, for many scattering processes 
within the dispersive band between states near the van 
Hove points, the bare interactions vanish, and the inter- 
actions away from these points and away from the Fermi 
level (which are kept in approximate form in our refined 
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FIG. 6: fRG results for a 32 x 3 discretization of the Brillouin 
zone, for the half-filled band, T = O.OOlt, V = O.Qt, J = OAt 
and U = 0. Left: Effective interactions VX(k, k', k + Q) with 
Q = (ir, 7r) and k near the (n, 0)-point at angle (solid line) 
or near the (n/2, 7r/2)-point at angle 7r/4 (dashed line) vs. k' 
moving around the FS at scale A = 0.33£ where the strongest 
attractive component of the coupling function exceeds — 24t. 
Right: Flow of the zero-frequency s-CDW (solid line) and 
d-CDW susceptibility for the same parameters. The inset 
shows the 32 FS points (filled circles) and the additional patch 
centers at band energy ±0.4f (open circles). 

patching scheme) are apparently not sufficient to produce 
an instability at a higher scale. 

For nonzero nearest-neighbor interactions V and J the 
fRG produces run- away flows at measurable scales. For 
instance, for U = 0, V = 0.6?; and J = OAt we find a 
flow to strong coupling, i.e., a divergence of some com- 
ponents of V/\.(k, k', k + q) in the one-loop flow, at scales 
A c « 0.033i. This scale sets an upper temperature bound 
for possible long range order. The divergence is dom- 
inated by scattering processes V\(k, k',k + Q) with k 
and k' at the FS and with wavevector transfer Q = (it, it) 
which flow to strongly negative values. For incoming and 
outgoing particles near the saddle points, these processes 
belong to the gi-type scattering known to diverge to — oo 
from the previous section. In the two-patch model, also 
the 52-type processes seemed to diverge. This tendency is 
not seen in the fRG for the full Fermi surface, now we ob- 
serve only a small growth of the <?2-type processes. In Fig. 
[6] we display the diverging component Vk(k, k',k + Q) 
of the effective interactions for two different choices of 
k vs. k' moving around the FS at a low scale A where 
the maximal component of the coupling function exceeds 
—2At. One clearly observes that the data for k fixed near 
(ir, 0) is basically a superposition of an attractive offset 
and a cos k\ — cos k^ modulation. The attractive offset be- 
comes more negative with decreasing scale and drives the 
s-CDW susceptibility. The d-wave part also grows and 
feeds the rJ-CDW channel. For k fixed near the BZ diago- 
nal, i.e., near the nodes of the cos k\ — cos ki form factor, 
the d-wave part is missing or strongly deformed, while 
the attractive offset is more or less unchanged. Hence 
we can expect a stronger growth in the s-CDW chan- 
nel. This is also clearly reflected in the data for the 
flow of the susceptibilities shown in the right panel of 
Fig. [5] The growth of the s-CDW susceptibility by far 
exceeds the d-CDW channel and all other channels such 
as spin density waves. Hence the fRG treatment agrees 
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with the meanfield picture developed in the previous sec- 
tion. Both suggest a s-CDW ordered ground state in the 
effective one-band model, translating into a bond-order 
wave as primary order in the full checkerboard system. 

The s-CDW divergence described here occurs in a 
rather large parameter window. We increased the value 
of U from up to 2t and found an increase of the critical 
scale from 0.033£ to 0.037i together with an even stronger 
dominance of the s-CDW over the d-CDW channel. 
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APPENDIX A: CHECKERBOARD LATTICE 



IV. CONCLUSIONS 

We studied the Hubbard model including n.n. repul- 
sion and n.n. antiferromagnetic exchange for strongly 
correlated electrons at quarter-filling (n = 1/2) on the 
checkerboard lattice in the strong-coupling (U = oo) 
and in the weak-coupling regime. The MF calculations, 
the ED simulations and the RG analysis all predict a 
spontaneous breaking of the translation symmetry in the 
ground state, that is characterized by the nesting vec- 
tor Q = (ir,n). The ground state is twofold degenerate 
as there are two types of crossed plaquettes ("strong" 
and "weak" plaquettes). These ground states can be 
approximatively described as product states of localized 
two-particle states on the "strong" plaquettes, which are 
equal amplitude superpositions of a singlet wave function 
on each of the six bonds. In this way the two electrons 
forming a singlet gain the magnetic energy — J and the 
kinetic energy — At. Furthermore, the presence of a third 
electron on the crossed plaquette which would cost a n.n. 
repulsion energy 2V is avoided. Therefore, the BOW in- 
stability is a cooperative effect that results from a simul- 
taneous optimization of kinetic and interaction energy. 
Note, that for negligible kinetic energy (V ^> t) other 
phases can be realized on the quarter-filled checkerboard 
latticed An interesting issue for further theoretical stud- 
ies are the properties and the excitations of the weakly 
doped BOW phase, e.g., the question whether the ad- 
ditional holes or particles pair and lead to a supercon- 
ducting phase where the U(l) gauge symmetry and the 
translation symmetry are simultaneously broken. 

In conclusion with the discovery and the description of 
the BOW phase on the quarter-filled checkerboard lat- 
tice we show how a strongly frustrated system can avoid 
the local frustration at fractional filling by exploiting the 
bipartite arrangement of larger units. We hope that with 
our work we can stimulate further progress in the chal- 
lenging and exciting field of strongly correlated electrons 
on frustrated lattices. 
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The elementary unit cell of the checkerboard lattice 
contains two lattice sites situated at (y = 1,2) 

where the vectors x y are two primitive lattice vectors. 
With this convention we choose the origin of the lattice 
at the center of a crossed plaquette, where there is no lat- 
tice site. The tight-binding Hamiltonian for this system 
is given by 



H„ = -t 



V M (c 



(Al) 



c r-xi+x„,i/cr) + n ' C ' 



where v = 2, 1 if v = 1,2, N is the number operator 
and fi = — 2t in the following. We introduce Fourier 
transformed operators as 



(A2) 



where N is the number of unit cells. With these operators 
the tight-binding Hamiltonian reads 



(A3) 



where k v = k • x„. Diagonalizing this Hamiltonian leads 
to a flat band at At and to a band with the dispersion 
£k = — 2t cos fc„ which is nothing but the nearest- 
neighbor tight-binding dispersion of the square lattice. 
The operators of this dispersive band are expressed in 
terms of the original operators by 



7ko 



(A4) 



with rk = 5Z„ cos 2 (£v/2). In weak coupling, we can re- 
strict our attention to the states close to the Fermi sur- 
face, where rk = 1 and for every operator on the checker- 
board lattice we can obtain an effective operator on the 
square lattice by the substitution — ► cos(/c I/ /2)7k CT - 
In this way the effective Hamiltonian © is obtained from 
the interaction Hamiltonian 
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